The Structure of a Vortex in the t-J Model 
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We study the single- vortex solution of the t-J model within resonating- valence-bond (RVB) mean- 
field theory. We find two types of vortex cores, insulating and metallic, depending on the parameters 
of the model. The pairing order parameter near both cores have d x 2_ y 2 +ir/d xy symmetry. For some 
range of t/J the calculated tunneling spectrum of the metallic vortex core agrees qualitatively with 
the STM tunneling data for BSCCO. 



Recent STM experiments JlLB] on vortex states in high 
T c cuprates have revealed many interesting properties 
that are not adequately described on theoretical grounds. 
Most of the existing literature relies on conventional 
Bogoliubov-de Gennes (BdG) aproach with a d-wave or- 
der parameter In the mean time it is widely ap- 
preciated that the high temperature superconductors are 
doped Mott insulators, and there is no a priori reason for 
the conventional BdG description to be applicable there. 

In the present paper we report the first unrestricted 
(see below) mean-field theory study for the vortex states 
in the t-J model with Coulomb interaction. The main 
results are as follows. 1) Depending on parameters of 
the model there exists two types of vortices, one with 
insulating core and the other with metallic core. At a 
fixed t/J and Coulomb interaction strength, the insulat- 
ing core is favored by low doping while the metallic core 
is favored by high doping. 2) Near the core of the vortex 
the pairing order parameter has d x i_ y i +if]d X y symmetry 
(see later discussions). The value of r\ tends to increase 
with doping. 3) The total integrated single-electron spec- 
tral weight is proportional to the local concentration of 
holes. (In reality the lost spectral weight will appear at 
energies above the charge gap which is taken to be infi- 
nite in the t-J model.) 4) The details of the tunneling 
spectrum inside a metallic vortex core, such as the exis- 
tence of zero-bias peak depends sensitively on the 
parameter choices for t/J, Coulomb strength, and dop- 
ing. However the gross feature that the coherence peak 
tends to be suppressed in exchange for low-lying spectral 
weight is observed in all cases. For the insulating core, 
the spectral weight is zero due to the Mott constraint. 
For both cores the background d-wave behavior is recov- 
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ered within a few lattice spacings from the center of the 
core. 

In standard notation, the Hamiltonian of the t-J model 
is given by 



(ij) 



In addition, due to the strong on-site repulsion the low- 
energy Hilbert space is constrained to have no more than 
one electron per site. 

So far this Hamiltonian has evaded exact solution in 
space dimensions greater than one. Limited exact di- 
agonalization results are often too small in system size 
to make a statement about complex experimental situ- 
ations. Meanwhile quite a lot is known about the var- 
ious mean-field states of this model (?|]|]. In particular 
the qualitative prediction of the phase diagram by the 
"RVB" mean-field theory || agrees with the experimen- 
tal findings. Recently one of us showed that in the su- 
perconducting state the gross prediction of the mean-field 
theory survives the low-energy gauge fluctuations || . 

In this paper we extend the mean-field treatment of 
Ref. H to study a single superconducting vortex. Since 
the presence of a vortex breaks the translational symme- 
try, we allow all the mean-field order parameters to be 
site/bond-dependent. We call such calculation an "un- 
restricted" mean-field theory. A similar approach has 
been used recently to study the stripe phase of the t-J 
model [[ll] . We believe that this mean-field theory should 
be adequate to describe the relatively high-energy and 
short-distance physics of the vortex core. 

Our starting point is the boson-fermion Lagrangian 
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(KijHij + h.c.) — mnj}. 
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In the above P tj = e aa i f irT fj a i , = f ja fi<j (summed 
over spin), Xi is the Lagrange multiplier that ensures the 
occupancy constraint, and Ay and Kij are Hubbard- 
Stratonovich fields. In the mean-field approximation 



bi, Xi, Kij, Aij are all treated as time- independent clas- 
sical fields that minimize the action. Note that we have 
also included the long-ranged Coulomb interaction in 
the model. Since we are treating b^s as classical fields, 
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the Coulomb interaction is in effect incorporated in the 
Hartree approximation. Making use of the invariance of 
the action under local gauge transformation bj — > e l ^ j bj 
and fj — > e^jfj, we will restrict bi to be real and pos- 
itive. The values of t/J and V c /J appropriate for the 
cuprate superconductors are not known exactly. In this 
paper we study a range of values for such parameters and 
several doping concentrations. 

The homogeneous mean- field ground state of the above 
action is the e?-wave superconducting phase first worked 
out by Kotliar and Liu . In our notation, such ground 
state corresponds to A^ = A, and bi = %Jx at every site 
(x, =hole doping). To test our calculation we first study 
the homogeneous system using the unrestricted mean- 
field theory. We find that for doping concentration x 
that is neither too small nor too large the solution is the 
uniform d-wave superconducting phase As a typical 
example, the electron density in the central 12x12 region 
is shown for t/J = 1.25, V c /J = 1.25 and x = 30/256 in 
figure 1(a). Due to the open boundary condition used, 
the electrons tend to accumulate at the edge (not shown) , 
leaving a somewhat higher hole density in the interior. 
We note that the electron density actually drops before 
it increases. When the same unrestricted mean-field the- 
ory is applied to the insulating limit, x = 0, we find 
the plaquette-valence-bond state first discussed by Af- 
fleck and Marston 0. The fact that we do not get the 
Neel state for zero doping is a well-known artifact of the 
RVB mean-field theory. At large enough doping (x > x c ) 
we find Aij — and the solution describes a Fermi liquid. 

We now turn to the case of a hc/2e vortex in the 
superconducting state. In the following we shall only 
present our 16 x 16 results for t/J — 1.25. We single out 
t/J = 1.25 in our study because in its neighborhood the 
calculated tunneling spectrum agrees qualitatively with 



the experimental findings. For considerably larger t/J 
the calculated tunneling spectrum deviates from the ex- 
perimental measured one. For example for t/J — 3 and 
x = 30/256 we find a zero-bias peak in the tunneling 
spectrum. It is clear that we are more interested in re- 
sults that are independent of the parameter choice. The 
rest of the parameters we use are 0.125 < V c /J < 1.25, 
and 16/256 < x < 38/256. 

The vorticity is imposed via the initial parameters 

Kij = K , = A° e ie - , b, = V£, A, = Mo (3) 

where -Ko,A?, /io are the bulk mean-field parameters 
(A? has d x 2_ y 2 symmetry), and 6ij measures the an- 
gle made by the position vector of the center of the (ij) 
bond and a fixed axis. We subsequently update the 
above parameters self-consistently. Upon reaching self- 
consistency we find that A,, = A^ e %&ii , with Ay having 
d x 2_ y 2 + irjd X y symmetry [ p0| . 

One important feature of our mean-field solution is 
that it obeys 

(bibi) + (f ia f i(7 ) = 1, Vz 

(Kij+tibjhMfivfjv) + tKyfabj) - c.c. = 0, W(ij). (4) 

These equations imply that, on average, the total boson 
and fermion 3-currents do not fluctuate. This is a direct 
consequence of the strong correlation physics inherent 
in the t-J model. In the field-theoretic treatment ||, a 
space-time local version of the above constraints arises as 
a consequence of integrating out the gauge fluctuation. 

A main conclusion of our work is the existence of two 
types of vortex cores - one with an insulating and the 
other with a metallic core. Small V c and low doping den- 
sity favors the insulating core, and otherwise the metallic 
core is favored. 
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FIG. 1. Electron density profile for (a) no vortex (x — 30/256, V c /J ~ 1-25), (b) vortex with metallic core (x — 30/256, 
Vc/J = 1.25), and (c) vortex with insulating core (a; = 16/256, V c /J = 0.125). Only the central 12x12 plaquettes of the 16 x 16 
lattice is shown. 



In Table || we summarize our findings for a number of 
different parameter choices. Based on this study we pro- 
pose that at a fixed Coulomb interaction strength, the 
vortex core changes its nature from metallic to insulating 



as the system is progressively underdoped. 

Vortex with Metallic Core: In figure 1(b) we show 
a typical example of a metallic-core vortex. We plot the 
electron density profile near the center of a vortex for 
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V c /J — 1.25 and x — 30/256. First we note a slight 
increase of electron density in the vortex core. In figure 
2(b) we show the local density of states (DOS) associated 
with a site on the central plaquette of the same vortex. 
Compared with the bulk DOS (black curve) there seems 
to be more states at low energies in the vortex core (red 
curve). As we step away from the central plaquette the 
two curves become indistinguishable beyond 4-5 lattice 
spacings. As mentioned earlier the self-consistent order 
parameter A^- near the vortex core shows d x 2_ y 2 + ir]d xy 
symmetry 0. For x = 30/256 and 38/256, r\ is about 
30% in the immediate vicinity of the core center. For 
x = 24/256 and x = 16/256, we find r\ « 20% and 3% 
respectively. Presumably a larger circulating current as- 
sociated with higher doping is responsible for this change 
in rj. The appearance of the id xy component inside the 
vortex and its consequence on the tunneling spectra were 
considered by several authors 

Inside the metallic vortex core, the low-lying DOS de- 



pends sensitively on the doping level. This is illustrated 
in figure 2 where we fix V c / J at 1.25 and vary the dop- 
ing level among x = 16/256,30/256, and 38/256. For 
x = 16/256, some of the spectral weights for \E\ > J is 
lost, while the lower energy portion is nearly unchanged. 
The integrated DOS is clearly reduced by the vortex. For 
x = 30/256, much of the lost spectral weight under the 
peaks at E = ± J has re-emerged at lower energies. For 
x = 38/256, the peaks at E — ±J are entirely gone, 
and the DOS near E = is considerably higher. We be- 
lieve that the low-lying DOS profile is determined by the 
following two factors. l)The presence of the id xy compo- 
nent opens up a gap and hence pushes the states away 
from E = 0. 2)The circulating current Doppler-shifts (ll| 
the quasiparticle energy levels and increases the DOS at 
E = 0. The dependence of the detailed shape of DOS on 
the choice of t/J,V c /J,x simply reflects the variation of 
the above two factors with the parameters. 
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FIG. 2. Density of states (red) at the center of the metallic vortex core for x = 16/256, 30/256, and 38/256 holes and 
Vc/J = 1.25. Dark curves represent the pure d-wave results for the same parameters. 



When comparing the local DOS in the vortex core 
and far away, it is important to bear in mind that in 
the present description the total integrated electron local 
DOS is proportional to \bi\ 2 . Since the vortex core has a 
higher electron density compared with the bulk, the cor- 
responding \bi\ 2 is smaller. As a result it will appear that 
some electron spectral weight has simply disappeared in 
the vortex core. In reality, the lost spectral weight should 
appear at energies greater than the Hubbard gap, which 
is treated as infinity in the present model. 

The extra low-energy DOS induced by the vortex di- 
minishes with distance from the core. Experimentally 
Pan et al. have measured the excess conductance at a 
fixed bias voltage in BSCCO and found that it decays 
exponentially For the system with x — 30/256 and 
V c /J — 1.25 (Fig. 2(b)), we have calculated the extra 
DOS at a fixed energy, E = 0.3J. Figure ||(a) shows 
such "excess conductance" for the center 10 x 10 sites. 
We have also plotted the conductance along the diagonal 
(nodal) direction with distance in Fig. 3(b), where the 
straight line is a pure exponential behavior. The fall-off 
distance deduced from this plot is approximately three 



lattice spacings. The excess conductance shows consid- 
erable amount of angular variation, with the maximum 
occurring along the nodal directions. Such anisotropy 
was not seen within the experimental resolution Q . 

Clearly we do not intend to compare our results quanti- 
tatively with the measured ones. It is however significant 
that there is a range of (reasonable) parameter choice 
that yields results in qualitative agreement with the ex- 
perimental findings. For this reason we believe that the 
model we study does capture the essence of the vortex 
physics in the cuprates. 
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FIG. 3. (a) Tunneling conductance difference between vor- 
tex and non-vortex states. The central 10 x 10 lattice is 
shown where each site is represented by a plaquette (dual lat- 
tice). Bright (dark) areas have a larger(smaller) conductance 
enhancement for the vortex, (b) Semi-log plot of the conduc- 
tance difference versus distance (lattice spacing=l) along the 
nodal direction. The linear fit is made with a slope of —0.34. 
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Vortex with Insulating Core: In figure 1(c) we 
show a typical example of an insulating-core vortex, 
found for V c /J — 0.125 and x — 16/256. In this case 
the electron density in the central 4x4 plaquette nearly 
reached one. The coupling of the central insulating re- 
gion and the surrounding is very weak, i.e. the links 
connecting any of the core sites with the non-core sites 
have Ay w 0, Kij « 0. Since \bi\ 2 w in the core there is 
negligible single-electron spectral weight there. The spin 
excitation spectrum in the vortex core is almost the same 
as that of the plaquette-valence-bond state of the insu- 
lating limit, i.e. a spin gap approximately equal to 2J 
separates the ground state singlet from the first excited 
state. It is possible that this feature is an artifact of the 
inadequacy of the RVB mean-field theory to correctly de- 
scribe the N(e)el state. Outside the vortex core \bi\ 2 > 
and the electron local DOS gradually recovers its bulk 
value. As in the metallic case the integrated spectral 
weight suffers a factor of \bi\ 2 /x suppression compared 
with the bulk. In the vicinity of the core we observe 
d x 2_ y 2 + ir)d xy order parameter. For underdoped materi- 
als, the evolution of DOS in the above-mentioned manner 
as the tunneling tip moves away from the center of the 
vortex will be an indication of the existence of insulating 
core. The possibility that for low doping the vortices can 
have insulating cores has been discussed in several earlier 
works ||Jl2|,[l6| . The existence of this type of vortex will 
be a clear manifestation of the proximity of the supercon- 
ducting system to the Mott insulator. We believe that 
the insulating vortex core and the appearance of charge 
stripes in the underdoped cuprates have a common ori- 
gin - at a particular doping density the system is on the 
brim of phase separation ||l7| , pT| . 

We conclude by noting that, since low-energy- 
integrated spectral weight is directly proportional to the 
local hole density, a careful measurement of its varia- 
tion can, in principle, reveal the local charge distribution. 
This suggests an interesting possibility of determining the 
presence of a vortex or a charge stripe by direct imaging 
in STM experiments. 
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TABLE I. Nature of vortex cores for various dop- 
ing (horizontal entry) and Coulomb interaction strengths 
(I=insulating, M=metallic). 
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